-it 


•v 

V * 



« I 


X-591-73-351 

PREPRINT 


urn tu x-. >^ss-r 


CORRECTION OF LASER RANGE 
TRACKING DATA FOR ATMOSPHERIC 
REFRACTION AT ELEVATIONS 


ABOVE 10 DEGREES 


j. W. MARINI 
C. W. MURRAY, JR. 


NOVEMBER 1973 


» W«-» 

* H M** 

* M* 

U) CO Ci u> 

rs"? 

H HI W 
H W CM 
O' O •>* 9 

c,, aY 

T 1 *• M 

H a o 
SC Ci 

o w UJ 

4* c* 

*gi < H 

• ^ ► 

O *4 £ 

CH?0 

o o w 

* W 5 W 
tn r»i 
r n 

133 3H 

o oo 
< 10 * 
W *3 
SB o 

o sr 

r. nr 
w o o ► 
o w 52 

t* CJ w 
to w 



GODDARD SPACE FLIGHT CENTER 

OREENBF.lt, MARYLAND 



N> a 
HO 
o o 


as P> 

to 


as 

H 

X* 

I 

-* 

—A 

tn 


CORRECTION OF LASER RANGE TRACKING DATA FOR ATMOSPHERIC 
REFRACTION AT ELEVATIONS ABOVE 10 DEGREES 


J. W. Marini 
C. W. Murray, Jr. 


November 1973 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 



» > 

<• 


CORRECTION OF LASER RANGE TRACKING DATA A ™ OSPHERI ° 
REFRACTION AT ELEVATIONS ABOVE 10 DEGREES 


j. W. Marini 
C. W. Murray, Jr. 


ABSTRACT 

A formula for correcting laser measurements of satellite 
range for the effect of atmospheric refraction is given. The 
corrections apply above 10* elevation to satellites whose 
heights exceed 70 km. The meteorological measurements 
required are the temperature, pressure, and relative humi- 
dify of the air at the laser site at the time of satellite pass. 

The accuracy of the formulr was tested by comparison with 
corrections obtained by ray-uacing radiosonde P™ files ; 
The standard deviation of the difference between the refrac- 
tive retardation given by the formula and that calculated by 
ray-tracing was less than about 0.04% of the retardation or 
about 0.5 cm at 10° elevation, decreasing to 0.04 cm near 
zenith. 
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CORRECTION of laser range tracking data for atmospheric 
correctioN r OF^las at elevations above 10 dkckebs 


INTRODUCTION 

The correction of tracking data for atmospheric refraction has been c^mively 
studied and many correction formulas have been published 1161. tor cer 
Sl^ceTn Physics applications, however, position accuracies of better 
thaVa few centimeLrs are desirable 17] , and these accuracies are much greater 
tan required for most previous applications. Out of the work cited, only ta 
approach given by Marini 131 . and the expansion and integral evaluations of 
Saastamoinen 15,01 provide the desired accuracy at lower elevation angles (10P 
Kfs. report Saastamoinen's integral evaluations are incorporated into 
Marini’s continued fraction form to provide relatively simple algorithms 
Erecting laser range-data using surface meteorological measurements. 


( 2 ) 


REFRACTIV1TY AT OPTICAL FREQUENCIES 

There are a number of formulas (8-111 for the refractive index n of air and for 
the corresponding refractivity 

N s 10 6 (n * 1) W 

all of which have sufficient accuracy for use here. The formula employed is (121 

/ 1.6288 0.01 36\ / P \ ( ! .) 

N = (287.604+ — + ^4 J y J 0 1 3.25 / \1 +0.003661 t / 

'° 055 (lOTOr) (r+ 0.003667) 

where 

X = wavelength of radiation in microns 
P = atmospheric pressure in millibars 
e ■ partial water vapor pressure in millibars 
t = temperature in degrees Celsius 

Because air is dispersive at optical frequencies, the group refractivity N g is 
also required 

d .... .. . dN 

Vdf<«>- N -*3r 


i 


(3) 


where f is the frequency. The expression for the group refractivity can be 
written as 

N g * 80.343 f(X>Y * 1».3 y 


where 

P = Total air pressure in millibars 
e = Partial pressure of water vapor (mb) 
T = Temperature (°K) 

and 


f(A) 


s 0.9650 + 


0.0164 0.000228 

X 2 X 4 


which, at the 0.6943 micron wavelength of the ruby laser becomes 

f(0.6943) = 1.0000 


(5) 

( 6 ) 


GEOMETRY AND NOTATION 

The geometry of the satellite-tracking station configuration is shown in Figure 1. 
Spherical symmetry is assumed, i.e. the refractivity is taken o e a _ 
height only. The height h is measured from the tracking station upward. Th 
subscript "0" designates quantities evaluated at the tracking station, the subscript 
"1", quantities evaluated at the satellite. The ray or phase path between track it^ 
station and satellite is shown as a curved line. 1.- Tue range R to the dis 
along the straight line connecting the tracking station and the satellite, and the 
tr^/elevation angle E is the angle between this line mid the h ^ z °^ la ^ th . t 
station. The nominal earth radius used is r e =• 6378km, and H is the height 
ofthe tracking station above sea level. The latitude of the tracking station is* 

degrees above the equator. 


EXPANSION FORMULA 


The apparent range R e 
by a pulsed system is 
the phase path (13,141 


between the ground station and the satellite as measured 
given by the integral of the group index of refraction along 


2 



(7) 


r "» 

*• * / sin# 

•''ft 


dr 


where the angle # Is given by Snell's law lor a spherically stratified medium 

nr cos# * n 0 r 0 cos# 0 <8) 

The oorreotlon sought is the difference between the measured and the true value 
of the range 

AR = 1% - R (9) 

The expansion ol HR in Inverse powers ol sin# 0 . following Marini 131 gives 

ARi 10* 

j*hN g dh-IO* 12 N 0 j N # dh 


( 10 ) 


♦ lOr^JlN^-l^ldh]-^ 


+ 000 

where ihe range '^rsr^rsatm^r^t 

Sit SEST«* - 1 neglected? as sbownln Appendix 1. because «10> 
is to be applied only where E > ltf and h, > 70 km. 

The — (10) is not the most useful one for nurny orbit determination pro- 

expansion of the angular correction is used 

0 O * E - 1<T 6 N 0 cot E (ll) 

substituting (11) into (10), and making suitable approximations 


4 


* 


AR= [|0'‘/VU,]--L 

‘ [v /"V* l 0 r,! I < N vi N2)Uh ] 

+ 000 

Equation (12) above is the expansion that provides the basis for the correction 
formula that is the subject of this report. 


EVALUATION OF INTEGRALS 

The evaluation of the integrals, appearing in (12), as functions of the pressure, 
temperature, and relative humidily of the surface air at the tracking station, 
has been treated by Saastamoinen 16] . For completeness, and because they 
differ in detail, our evaluations are given in Appendix 2. The results are 

KT 6 f N g dh * tO.O02357P 0 + 0.000141e 0 ] (13) 


~ I hN dh * f(X) (1.084 x ltf 8 ) P 0 T 0 K (14) 

r 0 J ® 

p 2 2 

l(T i2 J (NN g -yN 2 )dh = «X) (4.734x1 Or 8 )~ • (15) 

where 

f(*,H) = 1 -0.0026 cos2*>- 0.0003 1H (16) 


and 


K - 1.163 - 0.00968 cos2tp 
- 0.00104 T 0 + 0.00001435P 0 


(17) 


i 

i 


\ 


1 


4 

J 


\ 




CORRECTION FORMULA 

The formula for calculating the range error AR from the satellite elevation h 
is obtained by approximating (12) hy a continued fraction form 


AH = J (X) - 

AR f(*.H> 


sinli t 


At B 

" B/( A t B) 
sin I'' 1 0.01 


(IX) 


{ 

1 

i 

i 

i 


I 

\ 


where 


A = 0.002357P,, t 0.00014 lc 0 

P 0 2 2 

B s (1.0X4x 10’ 8 )P 0 T„ K + (4.734 x 10 H ) ^ - l/K) 

K = 1 . 1 63 - 0.00968 cos2<p - 0.00 1 04 T () + 0.0000 1 435P 0 


(19) 

( 20 ) 

( 21 ) 


Here 


AR 
E = 
Po 85 
T 0 = 

e o = 
f(X) - 

f(<A H) = 


= Range correction (meters) 

= True elevation of satellite 

Atmospheric pressure at the laser site (milliba ) 

Atmospheric temperature at the laser site (degrees Kelvin) 

Water vapor pressure at the laser site (millibars) 
l for a ruby laser, and is given by (5) otherwise 

1 tor a laser site at latitude and at sea level, and is given by (10) 
for sites at different latitudes <P and elevations H (in km) 


The water vapor pressure e 0 may he calculated (rom a relative humidity mea- 
surement R h (%) 

7.5 (T 0 -273.15) 


e„ * 


_?!l 

inn 


x6.ll x 10 


237.3 + (T 0 * 273.15) 


( 22 ) 



: (18) because toe error thereby in- 

:r,n l xrxzszz - 

SS iTCiS VES a small bias that oeeur. to toe expansion (12, 
ecause of approximations made in its derivation. 


iST OF ACCURACY 


6 


obtained by ray-tracing radiosonde refractivity profiles. 1 ho ray-ti aco cor- 
rections are considered to have stato-of-the-art accuracy, so that the differences 
between these corrections and those calculated from the simpler formulas re- 
present the penalty paid for simplicity in calculation and measurement. 

The data used in Figures 2-11 was obtained from the National Climatic Center at 
Asheville, North Carolina. It consists of radiosonde observations taken near 
Dulles Airport, Virginia, during the year 1967. 

Using the procedure described in Appendices 3 and 4, 634 refractivity profiles 
were calculated up to a height of 1000 kilometers from the rad osonde observa- 
tions. The calculated profiles were ray-traced 1161 at arrival angles of 1CT , 

15° 20° , 40° , and 80° , and the tropospheric errors in range and elevation angle 

were obtained. The histograms of these errors are snown in Figures 2,4, 6, 8, 
and 10. The correction formula (18) was applied using only surface data and the 
known elevation angle to obtain approximate tropospheric corrections. The dif- 
ferences between these algorithm corrections and the ray-trace corrections were 
calculated. The histograms of these differences is rhown in Figures 3, 5, 7, 9, 
and 11. The maximum bias of the error remaining after correction was -0. l cm, 
and the maximum standard deviation was 0.49 cm at 10° , decreasing to 0.04 cm 

at 80° . 

In addition formula (18) was compared with range corrections obtained by ray 
tracing (at arrived angles of 10% 15% 20% 40% and 80') radiosonde refrac- 
tivity profiles calculated at Jananarive (85 profiles), Fairbanks, Alaska (200 
profiles), Athens, Georgia (200 profiles), Greensboro, North Carolina (200 
profiles), and Nashville, Tennessee (135 profiles). The maximum standard 
deviation of the error in the algorithm at 10° was 1 centimeter andthema^mvmi 
at 80° was 0.06 centimeters. The maximum mean error of the algorithm at 1(T 
was 0.16 cm and the maximum at 80° was 0.07 cm. 


CONCLUSIONS 

An equation that corrects laser range data for atmospheric refraction using sur- 
face meteorological measurements has been deiived, and a ‘' om ^f i80 “ ,™ a ® 
between the corrections calculated using this equation (equation 18) and the cor- 
rections calculated by ray-tracing through a radiosonde profile. The compari- 
son (Figures 2-11) indicates that the differences between the corrections calcu- 
lated be the two methods are negligible for practical applications. Hence 
accurate refraction correction of laser range date can be made witeout the re- 
quirement for radiosonde measurements or lengthy ray-tracing algorithms. 


7 
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Figure 8. Tropospheric Range Error at About 15 Degrees Elevation for Laser Frequency of 0. 6943 Microns 
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Tropospheric Range Error at About 10 Degrees 
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It ahould be notated out that only the relative accuracy of the two procedures has 
been tested and that errors caused bv factors common to both methods arc not 

W tor the group refractive Index la used 

both in (18) and In the ray-trace equations, and any error In Its nuteltude would 
reflect equally to the corrections. Similarly, the hydrostatic equation used in 
equation(2-l) and hence (18) is also implicit in the ray-tracing method because 
tbe heiffhts that annear to radiosonde profiles are not measured quantities but 
rather are calculated from the measured pressures, temperatures, and [ r ®J ati ^® 
humidities using the hydrostatic equation. Also, both methods assume horizo 
homogeneity. Saastamotaen t6l has estimated the standard error from such 
sources to be less than 1 or 2 centimeters at 10° elevation. 
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NEGLECT OF SATELLITE RANGE 
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APPENDIX 1 

NEGLECT OF SATELLITE RANGE 


The correction 


AR 


■f 


sinfl 


dr - R 


( 1 - 1 ) 


can be written as 


AR » 10‘ 6 


(V, 

f . 1 

] r sin0 
,r o 

J. sin0 d R 

L^o J 


( 1 - 2 ) 


The expansion of the first term in (1-2), using suitable approximations (31 , gives 

k 

j N g (N 0 - N) d! 


sin 3 0 t 

-l0r« 2 


+ o o o 


The expansion of the bracketed second term in (1-2), which represents the dif- 
ference AR_ between the geometrical lengths of the phase and the straight-line 
paths between the satellite and the tracking station, can be obtained by expanding 
equation (A5) of reference l3l in inverse powers of sin0 o giving 


AR, = — 5 10*' 2 / N 2 dh 

sin J 0 o 2 J 

1 (/ N dh) 2 eos 2 0 o 

— 10* 2 r 

2 R sin0 o 4 


(1-4) 


Al-1 


i 


1 


i 


\ 


\ * 
\ 





between the magnitudes of N amt Ng 

1 |0’ 6 /N dh cos 2 0„ n 5) 

relative error = - ~ ^ 

The satellite height h, is roughly approximated by approximated byR si n 0„. 
and the zenith integral is about 2 meters: 

relative error * tun 2 ( 1 ) 

. u ia «. e satellite height in meters. Taking hj ^ 70km and ^ 
error calculated from (1-6) is less than 0.05% which can be neglected. 



APPENDIX 2 


EVALUATION OF INTEGRALS 


From the perfect-gas law, the law of partial pressures, and the hydrostatic 
oquation 


dP _ Mg(P-e) Mwge 
' dh RT RT 

( 2 - 1 ) 

MgP 0.378 Mge 
= ' RT ' RT 

where 1151 

M = 28.966 = Molecular weight of dry air 
Mw = 18.016 = Molecular weight of water vapor 
R = 8314.36 Joules (“K)' 1 (Kg - Mole)' 1 
= Universal gas constant 
g = acceleration of gravity (m/s) 
h = height (m) 

Combining (2-1) and (4) 

f N g dh = - 80.343 f(\)-^-JjdP 

( 2 - 2 ) 

+ [30.5 f(X)- 11.3] J^rdh 

The first integral on the right side of (2-2) above can be evaluated using the ap- 
proximation ll 5 1 . 


g = 9.806 1 1 - 0.0026 cos2 <p - 0.00031 (H + h)] 

i. L- 1 1 + 0.0026 cos2i + 0.0003 1 (H + h)l 

g 9.806 1 


(2-4) 


•* rm 


M 


% • . 


i 


) 

t 



\ 

I 

1 


i 


from which, integrating the last term by parts, 



._L_- [ ) + 0.0026 cos 2<p + 0.0003 1 (H + jp“ / MMl 

= p 0 /i 


(2-5) 


where g is the value of g at the height 



above the tracking station or H + h above sea level. Saastamoinen uses a gravi 
iational constant evaluated at* 

H + h ■ 7 .? + 0.9 H km * 2 ' 7 * 


From (2-7) and (2-3) 

g - 9.784 ( 1 * 0.0026 cos 2<p - 0.00028 H) 


( 2 - 8 ) 


s 9.784 f(«p,H) 


where B is the station elevation in kilometers. Saastamoinen has also evaluated 
the integral 



where the g appearing in (2-9) is set equal to g in (2-8) as a convenient approx 
imation. The expression for the zenith integral becomes 


J Ndh - 80.343 
+ 130.5 «CX> - 1 1.31 4§f ®o 


*An equivalent result canbe obtained by numerically estimating h using (2-17) with T 0 set equal to 
T e + where T e is the sea level temperature. 


A2-2 


1 


10* 6 / N dh = 


f(X) 

K*,H) 


0.002357 P 0 


(30.5 • 1 1.3/f(X)) 

* \tt 


0.000141 c 0 


( 2 - 11 ) 


Neglecting small errors in the second term of (2-11), equation (13) results. 


SECOND INTEGRAL 

In equation (12), the magnitude of the coefficient of 1/sin E is ^^2.4 meters, 
while the coefficient of 1/sin 3 E, is only about \ centimeters. At E = 10 . the 
magnitude of the first term is about 12 meters, while the second is about half a 
meter. Consequently the second term need not be as accurately evaluated as the 
first, and it is sufficient to use the approximation 


Hi f hN dh = — f 80 343 f( * - hdh (2-12) 

r o J 8 r e J T 


where r e is a nominal earth radius (6378km) and the air is assumed to be dry. 
It is also sufficient to treat g as a constant throughout. 

From (2-1), and integrating by parts 

fi Mh = jr 8 / Mh <2 ' l3) 


The pressure P in (2-13) is obtained by integrating (2-1) 

p - p o e * p [^/? dh ] 

The temperature T is assumed to have a linear slope 


T * T 0 + 0h 

and the integration in (2-14) is carried out giving 


" Po 


yMg/R0 

; 


(2-14) 


(2-15) 


(2-16) 


A2-3 



'■"-f 

i. 


The integration in (2-13) may now be performed 


/ P dh = P 0 • Mg * i . jj£ 

Mg 


(2-17) 


mi 


From (2-12), (2-13), and (2-17) 


iSlThN, 

r o J 


.. • 10- 6 (80.343)R 2 p T * 

dh = f(X) xi r ° lo 


r e M 2 g 2 


= f(X) (1.084 x 10* 8 ) P 0 T 0 K 
where g has been set equal to 9.784 and the factor 


1-R0 

Mg 


(2-18) 


(2-19) 


is equal to unity In an Isothermal atmosphere (0- 0, and * » equri to abort »;»_ 

In an atmosphere In which the temperature lapse rate Is a constant 6 /to <P 

-6° Am). 

Rather than use the theoretical value for K given by (2-19), which is * 

« * lanao rate the value of K used in the corrections equations is taken 
hHn empirical constant which was determined by solving (2-18) for k tod to - 
culattngUsvalue by numerically Integrating through the atmosptoresof theUS. 
Stewart Atmosphere Supplements, 1966. Using linear regression on the values 

so obtained, the formula 

* K= 1.163 - 0.00968 cos2v> - 0.00104 T 0 + 0.00001435 P 0 (2 20) 

resulted. Here is the latitude of the tracking station. 


THIRD INTEGRAL 

The contribution from the third Integral In (12) Is only marginally significant, 
and the term can be approximated by 

ll0 « 2 f N g N dh = —-10’ 12 (80.343) 2 f(X) dh (2 ‘ 2 

2 a. 


A2-4 


( 2 - 22 ) 



Assuming a constant temperature gradient, and using (2-16) 


io- 


12 


J (n Ng -y N-) 


dh il!^ (8 0.343>’«X>^ 


1 

* 1 + RP 

2mg 


The last factor in (2-22) can toe expressed in terms of K using (2-19), giving (15) 

\ 


AM 



APPENDIX 3 

PROGRAM FOR CALCULATING REFRACTIVITY PROFILES 
FROM RADIOSONDE DATA* 


RADIOSONDE DATA 

Radiosonde observations are measurements of pressure, temperature, and 
humidity taken from the surface up to the point where the balloon that carries 
the sensors bursts 111 . The values of temperature, pressure and relative 
humidity measured at certain standard and significant levels during each balloon 
ascent from numerous weather stations is available from the National Chmatic 
Center. This data can be used to construct continuous refractivity profiles from 
the surface up to the point of highest measurement. Above the latter point, the 
refractivity profile can be extended by assuming a suitable temperature profile. 

GEOPOTENTIAL ALTITUDE 

The equations used to calculate the refractivity profiles employ the geo- 
potential altitude H 12, p. 2171 , which is given by 


1 f 2 

H = G J 0 


gdZ 


( 1 ) 


where Z is the geometric altitude, and the lower limit of integration is from ^ sea 
level (Z = 0). H is in geopotential meters when G equals 9. 8 m/sec . The local 
acceleration cf gravity is calculated from the latitude tf»by {2, p. 488] 


g 0 = 9. 780356 (1 + 0.0052885 sin 2 ~ 5.9 * 10 -6 sm 2 24>) 
and 12, p. 2171 


go r o 


g 


(m/sec 2 ) 


( 2 ) 


(3) 


( r o + Z) 2 

Here r 0 is an effective earth radius given by 12, p. 2181 

2 g 0 


r ° 3.085462 x 10" 6 + 2.27 * 10" 9 cos 2<t> - 2 * 10“ 12 cos 4<*> 


(m) (4) 


•This appendix is self-contained. It has separate references, and the notation used differs from that in the 
rest of the report. 


A3-1 


• • • 

* r 

£• i- 


From (1) and (3) the conversion between geopotential and geometric altitude 
given by 1 2, p. 218) 

io / r 0 M 
H = 0 


J 


go r o 


VIRTUAL TEMPERATURE 

The calculations also make use of the virtual temperature T, (3) which Is 
related to the ordinary temperature T (°K) by 


1 - 0.379 p 


where e Is the partial pressure of the water vapor In the air, and is given by 
)4, p. 343) 


( p \ 7.5 ( T - 273.15) 

(6.11) 10 237 - 3+ < T ' 2731S) (mbar) 


R h being the relative humidity in percent.* 
CALCULATION OF GEOPOTENTIAL ALTITUDES 


*tf the dewpoint temperature T d OO is given instead of the relative humidity, e can be calculated 
from (8) by setting R h " 1CJ * nd T = T d* 


table of pressure, temperature, and virtual temperature versus geopotential alti- 
tude. The virtual temperatures at the given points are calculated from the 
measured values of P, T, and R h using (8) and (7). 

To calculate the geopotential altitudes, it is necessary to assume hydrostatic 
equilibrium [ 3] 


dP - -p gdZ 


(9) 


The density P is given with sufficient accuracy by [ 31 


MP 

9 " RT V 

The apparent molecular weight of dry air is taken to be 1 2, p. 289J 

M = 28.966 


( 10 ) 


( 11 ) 


and the universal gas.constant 1 2, p. 2891 


R = 8314.36 Joules (°K) _I (Kg-mole)“ l (12) 

Using the assumption that the virtual temperature is a linear function of geo- 
potential height between any two adjacent measured points H, and H 2 , (9) may be 
integrated with the use of (1) and (10) to give 


OM(H, - H|) 



(13) 


which may be written as 


H a 



X in ( 1 * 2 /P|) 
In (1 + X) 




(• 


X 

2 


+ 


X_ 

3 


• t 


(14) 
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where 


(15) 


<T VJ - T„,VT. 


vl 


is established. 


CALCULATION OF REFRACTIVITY PROFILES 

The radio refractivity N is given by the formula! [5, p. 7] 


N 



3.73 x 10 5 


e 

T 2 


(16) 


with Paid? expressed In millibars and T in degrees Kelvin. 

To calovdate N at a given height, i.e., to obtain a point at a 

profile, U is necessary ,o know the values of P, T and e at that height. These 

are obtained as follows: 

The height is converted t. a geopotential altitude by adding it to ^ 

vapor pressure e is calculated from (7). Substitution into (16) then gives the 
required refractivity. 

•The geopotential altitudes are computed at the self-consistency of the 

(R and M, for example) adopted. 


(K ana m* lor exwiipw; •wt'"'— 

t At optical frequencies (2) and (4) of the main text are used. 
N - 77.6 t k2 ) 

where the wavelength \ is in microns. 
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A listing of the FORTRANH program with a sample profile calculated from 
meteorological data taken at Dulles airport on 1 January 1967 is shown in 
Appendix 4. 

Also shown are the surface measurements of temperature, pressure and 
relative humidity, the tropospheric range error obtained from ray-trace (RANGE 
ERROR), the tropospheric elevation angle error, the tropospheric range error 
approximation (RANGE ERROR APPROX) obtained from using equation (18) of 
the main text, and the difference between the ray-trace and the approximation 
(RANGE DIFF) for arrival angles of 10°, 15°, 20°, 40° and 80°, 
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STATION * 93734 YEAR * 67 MONTH * 1 1-AY * * HOUR * 

PAMjP . 0.277020160 0* KILOMET RES __ _ 

-»«M «P WHH»<K2«VtTV ^ -^M*0«02«»-04- 
HF. IGMT » 0.049872630-01 KILCMlTERS 
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■•- HO^R - 

STATION * 93734 YEAR * 67 MCNTM * I CAY l 
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su» OF- REFRACT I-V ITY -* 0*a36W»03O0 64- 
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rnhfi SWf m ArWa » O W- 94 *144. 1 BARS 
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Pi.RCf.NT 
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